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Abstract 

The modification and amplification of the gluon angular distribution pro- 
duced along with hard jets in nuclear collisions is computed. We consider the 
limit of a thin quark-gluon plasma, where the number of rescatterings of the jet 
and gluons is small. The focus is on jet quenching associated with the forma- 
tion of highly off-shell partons in hard scattering events involving nuclei. The 
interference between the initial hard radiation amplitude, the multiple induced 
Gunion-Bertsch radiation amplitudes, and gluon rescattering amplitudes leads to 
an angular distribution that differs considerably from both the standard DGLAP 
evolution and from the classical limit parton cascading. The cases of a single and 
double rescattering are considered in detail, and a systematic method to compute 
all matrix elements for the general case is developed. A simple power law scal- 
ing of the angular distribution with increasing number of rescatterings is found 
and used for estimates of the fractional energy loss as a function of the plasma 
thickness. 
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1 Introduction 



One of the new predicted observables in nuclear collisions at RHIC energies (v^s ~ 200 
AGeV) is jet quenching [l]-[5]. While coUisional energy loss of a hard jet propagating 
through quark-gluon plasma has been estimated to be modest [1, 2]: dEcou/dx <^ 1 
GcV/fm, for accessible plasma temperatures, the induced radiative energy loss, as 
estimated [3, 4] using the Gunion-Bertsch [6] non-abelian generalization of the Bethe- 
Heitler formula for QED is expected to be significantly larger with dE^ad/dx > few 
GeV/fm. An uncertainty principle analysis in [7] leads to a simple estimate 

dE {kl) 

where (k^) is the mean transverse momentum of the radiated gluons. In Ref. [8] a 
model for non-abelian energy loss was developed and a class of multiple scattering 
diagrams lead to dE/dx ~ 1 GeV/fm independent of the path length. However, in 
BDMPS [9] gluon rescattering diagrams in the medium were shown to induce a nonlin- 
ear enhancement of the energy loss because {kj_) = fi'^L/X. The result for an incident 
very high energy jet penetrating a plasma of thickness, L, in which the mean free path, 
A = l/{ap), and the color electric fields are screened on a scale /x, is [9] 

— as{p_L) (xasii J . (2) 

Here {p\) is the jet broadening due to multiple elastic scatterings of the jet in matter, 
and a factor of order unity is suppressed. This estimate is derived for the situation 
when the radiation formation length is large compared to the target thickness and the 
target is thick compared to the mean free path: 

ALf = E/fx^ > L > A . 

In the other theoretical limit corresponding to Landau-Pomeranchuk-Migdal (LPM) 
[10] regime in QED, L ^ ALf ^ A, and the energy loss was shown to saturate as 
dE/dx oc as\fpE. 

The remarkable prediction that the non-abelian dE/dx increases linearly with the 
plasma thickness in nuclear collisions, suggests that the total jet energy loss could 
increase quadratically with the nuclear size. However, several practical problems limit 
the application of the above results directly to observable jet quenching in nuclear 
collisions. First, jets are produced in such reactions with high initial virtualities in 
hard pQCD processes, and thus induced radiation from secondary interactions must 
compete with the usual broad gluon shower accompanying a hard pQCD event. 

Second, the observable energy loss depends on the detailed angular distribution of 
the gluons. Unlike in electrodynamics where the energy of the final electron and pho- 
ton can be directly measured, the energy of a single parton is only measurable when it 
is kinematically well separated from other jets. This is because the observable hadron 
showers have wide angular distributions. Nearly coUinear partons in particular are not 
resolvable because they are absorbed into the coherent hadronic wavef unctions. There- 
fore, in order to test experimentally the non-abelian energy loss of jets, dNg/dyd'^k± 



2 



of the induced gluon radiation must be known. The angular integrated energy loss, 
dE/dx, is insufficient to predict its experimental consequences. 

Third, for finite nuclei L/\ < 10 is not very large, and a large fraction of the jets 
are produced in the nuclear "corona" with L/X < 3. Thus the form of the angular 
distribution for only a few collisions is needed. 

Fourth, a problem specific to nuclear collisions is that the very large transverse 
energy background due to minijets (the partons from which the dense quark-gluon 
plasma is made ) prevents the application of traditional transverse energy jet cone 
measurements of jets [11]. Therefore, jet quenching must be inferred from the system- 
atics of leading hadrons at moderate p_L ~ 3 — 10 GeV relative to well know spectra 
in elementary nucleon-nucleon reactions. This step requires the interface of the parton 
level shower calculations with hadronization phenomenology. 

This paper is part I of a study of jet quenching concentrating on the formal problem 
of computing the angular distribution gluon radiation from off-shell jets produced in 
a dense but thin medium. In the subsequent paper II [12], we encode the results of 
this paper into a nuclear collision event generator (HIJING1.36 [13] using the Lund 
JETSET [14] string fragmentation scheme for hadronization) and compute the hadron 
inclusive jet quenching pattern. Note that in a separate study [15], the transverse 
momentum dependence of the LPM effect in QED was discussed in another formalism 
and compared to the analytic results for QCD derived in this paper. 

Multiple interactions in QCD are comphcated by the color charge of the radiated 
gluons. This leads to interference between the production amplitudes and the rescat- 
tering amplitudes that makes it impossible to factor the cross section into a product 
of the elastic cross section and the gluon distribution. In BDMPS, a great technical 
simplification was achieved by considering only the angular integrated dNg/dcu, energy 

— * 

distribution. The freedom to shift the transverse momentum, k±, integrations was used 
often to reduce the multiple collision problem into tractable form. In addition, a large 

Ug = L/X ^ 1 and iVc ^ 1 approximations were employed to convert the problem into 
an elegant equivalent Schroedinger equation problem valid under the conditions that 

1/At < A < L < L/ = E/fj,'^ . 

An alternative path integral approach proposed in [16] obtained similar results (see 
also [15]). 

Unfortunately, the continuum approximations developed in [17] are not applica- 
ble to thin plasmas and more importantly they do not apply at all to the problem 
of computing the angular distributions. The main complication is that the angular 
distributions require a complete treatment of all gluon and jet final state interactions. 
In this paper, we simply follow the direct albeit tedious route of computing the sum of 
squares of all diagrams with finite number of rescatterings. In the thin plasma limit, 
the main simplicity arises from the fact that for < 3 the number of amplitudes 
2ns+i _ ]^ ig g^i^ gmall enough for practical calculations. 

Historically, the on-shell single rescattering case, Us — 1, was first considered by 
Gunion and Bertsch [6] (GB). They estimated the hadron rapidity distribution due to 
gluon bremsstrahlung associated with a single elastic valence quark scattering in the 
Low-Nussinov model. The important difference compared to the problem addressed in 
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this paper is that GB computed the non-abehan analog of the Bethe-Heitler formula 
while we are interested in the analog of radiative energy loss in sudden processes similar 
to the beta decay. The GB problem is to compute the soft radiation associated with 
a single scattering of an incident on-shell quark prepared in the remote past, i.e., 
to — — oo, relative to the collision time, ti. In our case, the radiation associated with 
a hard jet processes in nuclear reactions must take into account the fact that the jet 
parton rescatters within a short time , ti — to ~ Ra, after it is produced. The sudden 
appearance of the jet color dipole moment within a time interval ^ 1/E results in a 
broad angular distribution of gluons even with no final state scattering. The induced 
radiation cased by rescattering necessarily interferes strongly with this hard radiation 
as indicated in [18]. 

To model a thin plasma, we employ the static color-screened potential model as 
considered in Gyulassy-Wang (GW) [8]. The scattering potential has Fourier compo- 
nents 

Here Ti{Ta{R)T}j{FC)) — dahCB.DR/ Da is proportional to the second order Casimir 
of the Dr dimensional representation R of the jet parton, and Da = — 1- In 
principle, the i*'^ target could be in the Di dimensional representation with Casimir 
C2{i) so that Tr(Ta(i)Tfe(j)) = 6abSijC2{i)Di/DA- Note that in eq. (3) corresponds 
to Ai{q)/{-2iP°) in the notation of Ref. [8]. 

For a thermally equilibrated medium at temperature T, the color screening mass 
in pQCD is given by = AnasT"^. In addition, as in ordinary plasmas, no gluon 
modes propagate below the plasma frequency, cUpi ~ A*/\/3. In practice, lattice QCD 
calculations of /i indicate sizable nonperturbative corrections to the pQCD estimates 
for achievable temperatures. Therefore, we simply take here ~ a;p; ~ 0.5 GeV as a 
characteristic infrared scale of the medium. In perturbation theory, there is a relation 
between the screening scale /i and the mean free path A, 

iJ^/X Ri Ana^p (4) 

where p is the density of plasma partons weighed by appropriate color factors. Another 
important scale for our problem is the Bethe-Heitler frequency, 

(^BH = ^A*^A > p (5) 

in the dilute plasma approximation assumed here. 

This paper is organized as follows: in Sec. 2 we derive the spectrum of the gluons, 
radiated as a result of the sudden production of a high energy parton in a highly lo- 
calized wave-packet J{p) in free space. The large radiative energy loss, which we call 
"self-quenching", grows as AE ~ Eo\og{Eo/ p). The strong destructive interference of 
this initial amplitude with the subsequent induced radiative and cascade amplitudes 
is considered in Sec. 3. The main formal results of this paper are derived in that sec- 
tion. In Sec. 3.1 jet acoplanarity is reviewed for the case of multiple elastic scattering 
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without radiative energy loss. The formal structure of joint jet plus gluon probability 
density including multiple rescattcrings is derived in Sec 3.2. A systematic procedure 
for extracting all the relevant matrix elements in the eikonal approach is developed 
using a simple binary encoding scheme in Sec 3.3 to help automate the computation 
of phase factors and the color algebra. In Sec 3.4 the interference form factors and 
their dependence on the detailed plasma geometry is derived. These form factors con- 
trol how close the exact result is to the extreme factorization or the classical parton 
cascade limits. In Sec. 3.5, the physical interpretation of the relative conditional prob- 
abilities computed from the square of the sum of amplitudes and the wavefunction 
renormalization necessary to compute the gluon angular distribution as a function of 
Ug is presented. Sec. 4 illustrates in detail the single rescattering case. We rearrange 
the amplitudes into Gunion-Bertsch and gluon rescattering terms that reduce in special 
limits to the known cases. An additional form factor associated with the transverse 
profile of the target is derived that helps to identify the physical interpretation of the 
rearranged terms. In the simpler "broad" target case, those form factors drop out. 
The induced gluon bremsstrahlung (beyond self-quenching) is shown to be strongly 
suppressed relative to the classical limit due to the destructive LPM effect in QCD at 
small angles. However, the angular enhancement at large angles due to gluon rescat- 
tering remains a key feature. The more involved case of two scatterings is treated in 
Sec. 5. The analytical result shows a power law enhanced angular broadening of the 
gluon distribution relative to the case oi Ug — 1- This approximate power law scaling, 
oc (1 + Ca/Cr)^" , is used in Sec. 6 to extract an approximate independent reduced 
quantum modification factor of the gluon angular distribution. We apply this scaling to 
compute the angular integrated gluon energy spectrum, dl^""''' /du and the jet fractional 
energy loss AE^^^yE. The numerical results show surprisingly an approximate linear 
dependence of those quantities on the plasma thickness at least up to L = n^A ~ 3 fm. 
The jet energy loss remains small because the wavefunction renormalization factor also 
grows rapidly with Ug. As emphasized in the Summary these estimates are only lower 
bounds because multi-gluon emission is not treated in this paper. The generalization 
to multi-gluon effects is deferred to ref.[12]. Details of the Ug — 2,3 matrix elements 
and their rearrangement into the physical hard, Gunion-Bertsch, and gluon cascading 
components as well as the relevant color factors are recorded in Appendix A and C to 
ease the fiow of the presentation. In Appendix B we compare the Ug = 2 case to the 
approximation of BDMPS where some of the terms in the amplitude were neglected. 
In Appendix D we derive the curious color "wheel diagram" that arises in the Ug — 3 
case. 



2 Hard Radiation and Self-Quenching 

Consider first the amplitude for the production of a hard jet with momentum Pq 
localized initially near Xq = (to,0). We denote this color independent amplitude by 

Mj = J{P)e'^^'' . (6) 

The hard vertex is localized within a distance ~ l/Po, and the amplitude, J{P), is 
assumed to vary slowly with P on the infrared screening scale of the medium. This 
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spectrum of hard jets is of course softened due to gluon radiation to first order in Qg 
as shown in Fig. 1. 

Mj 



.to 



Mj X Mq 

.to 



k,c 



Fig. 1. The Mj hard jet production ampUtude and the Mj(g)Mo soft gluon radiation 
amphtude. 



The final state in that case is characterized by a gluon and a jet with momenta 

k] 

= (Po, P., Q± - - [(1 - z)P,\ ~ f^l' , - J , (7) 

where (■ ■ •) denote their 4-momentum and [• • ■] denote their light-cone momentum 
coordinates. 

Wc emphasize that the "z" direction in this paper refers to the unperturbed hard 
jet direction. In nuclear reaction, the jets are directed perpendicular to the beam axis 
in the frame in which its longitudinal rapidity vanishes. Therefore kx,ky, kz actually 
correspond to a 90 degree rotation about the y axis of the conventional frame in which 

— * 

the z axis is parallel to the beam. The initial "transverse" momentum Q± of the jet 
thus accounts for both its longitudinal rapidity as well as its azimuthal (symmetric) 
dependence about the beam. The energy loss, u, leads to a suppression of the high p± 
distribution of the fragments emerging from the jet. 

The gluon polarization in the A'^ — gauge is e{k) = [0, {2e± ■ k^)/ {zPq), el]. We 
consider here eikonal kinematics where 



u = Jkl + k\^ zEq 



2zEq ' 



/ - V5hZ-(1--)^o + ^3^. (8) 

The amphtude to emit a soft gluon and the jet neglecting spin effects can be ex- 
pressed as 

Mj®Mo = i [ d'P J(P)e^^-° ' • + 5\P -k-p)^ 

J (p ky 
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~ e'^°'°J{k+p)Mo{k,p) , (9) 
where A^o is the radiation amphtude that in the eikonal hmit is given by 

Mo = -2tg/-^-^e'''^'ic. (10) 

Here, c = Tc{R) is the SU {Nc) color generator Tc in the Dr dimensional representation 
of the jet parton with Tr{cc) — CrDr. 

A simple direct way to recover the phase factor and the l/k\ singularity of M.q is 
through the use time ordered perturbation theory as emphasized in [9]. This follows 
from integrating over the emission time of the gluon 

r dt e--l*l+^(%-^^)* = ^e^*°& , (11) 

where 

Ef - Ei ^ u; + po - Eo k\/2uj . 

This method is particularly useful when final state interactions are taken into account 
in the next sections. 

The probability distribution for the jet production is given by 

W3p(0) _ n^l 12 '^^P - /j^°Vd) ^ ri2l 

dU - (27r)32pO - ^ ^^\27r)32|p| ' ^^^^ 

where p'-^-' [p] is proportional to the invariant hard production probability density with 
no elastic rescatterings of the jet nor radiative energy loss. The pQCD jet distribution 
is approximately a power law l/p^^ relative to the beam axis and uniformly distributed 
in longitudinal rapidity. The power g ~ 3 depending on the beam energy and the p±_ 
range under consideration. However, in terms of our rotated variables with p^ ^ Px 

p^'\p)r^c/ipi+pir 

independent of our p^ that is parallel to the beam direction. We emphasize this point 
because unlike the usual energy loss problem of a beam passing through matter, the pe- 
culiar longitudinally expanding cylindrical geometry of the matter in nucleus collisions 
makes applications of simple formulas in the rotated frame somewhat unintuitive. 
The double inclusive jet plus gluon probability distribution is given by 

^6p(0) _ (0) . d'k d^p 

d l^rad - Prad{l^,P) (3^)32^ (27r)32pO ' ^^^^ 

where 

p^,{k,p)^TT\Mj®Mo\' (14) 

is the double inclusive probability density for producing a jet and a gluon in the case 
of no rescatterings. 
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As seen from (9,10) the double inclusive density factorizes for soft gluons into a 
production probability of an off-shell jet with k + p and a distribution of gluons. 

d^k TT^ kl- ' ^^^> 

where the dependence on p is only implicit through the kinematic boundaries. As 
in QED the radiation pattern associated with a hard jet corresponds to a uniform 
rapidity distribution along the jet axis with a broad transverse transverse momentum 
distribution that is uniform in log/c^. 

If the off-shell dependence of the source amphtude is weak and the gluon is soft in 
the sense 

-L\j^\^\+^^P+k)\^^^L^\J(\p + k\,p + k)\' , 
IP I \p + k\ 

then we can write 

p^^J,{k,p) « 2(27r)V°)(p + k)R^'\k) . (16) 

The gluon radiation associated with any hard pQCD process leads to what one can 
think of as "self-quenching" of the jet in the sense that the initially produced high 
energy jet parton loses a substantial fraction of its energy 

A^o = Cr^ — -j^zEo^ Cr—Eo log Eo//. . 17 

TT Jo Z Jjj? kj_ TT 

This occurs even in free space where there are no final state interactions. In practice, 
the large logarithm implies that multiple gluon emission must also be considered. This 
leads to a Sudakov form factor for the jet and a probabilistic Alterelli-Parisi parton 
shower [19]. The final state multigluon shower can be calculated most readily by one 
of the many Monte-Carlo event generators such as PYTHIA [20]. Those generators 
encode empirical "parton to hadron" fragmentation functions and thus allow the de- 
tailed study of the effect of parton showering on the final hadron distributions. The 
magnitude of the self-quenching effect, i.e., initial and final state radiation, is found to 
be significant and is essential to describe accurately the high p± data as shown in [21]. 
Our main interest here is to derive the extra quenching associated with gluon radiation 
induced by final state interactions of the jet in a thin plasma. 

3 Gluon Radiation with Rescatterings 

In this section, the self- quenching case is generalized to take into account possible 
multiple scattering of the jet and gluon following the hard jet production. 

3.1 Multiple Elastic Scattering 

The amplitude that jet undergoes Us elastic scatterings with momentum transfers qi — 
(0, Qiz, qi±) is given by [8] 

d^Qi 



= / d'P J(P)e^^^o / n (j^, (2PX(i)A(P + Q,^,)V(q^ 
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x5\P-p + j:q,){TM---Ta„AR)) . (18) 

where A(p) = (p^ — + ie)~^ is the jet propagator neglecting spin effects, and the 
cumulative momentum transfer Qi = J2]=i Qj- The invariant jet probability distribution 
together with Ug unobserved recoil partons in the target and no gluon radiation is to 
lowest order (without Sudakov factor) given by 

a -^\.Y\Mj I ;j(27r)32p0 (27r)32pO ■ ^^^^ 

The ensemble average over the location of the scattering centers is taken as in [8] 

(•••).— / ft {dx,^T{x,^) ^h^^e-'^eiz, - ^,_i)) (•••)■ (20) 

In a static medium, the average longitudinal separation between the centers is spec- 
ified by the mean free path, A = Aj. An expanding medium, on the other hand, 
has Aj > Aj_i can only be approximately treated with a constant mean free path. 
The transverse coordinates are distributed according to a normalized overlap function, 
T{x±) ~ l/(7ri25_)- We assume that the transverse size R± of the target (relative to 
the jet axis) is large compared to the range of the interactions, 1/^. 

In the high energy [p'^ = E ^ uj ^ fi) eikonal limit, going back to the picture 
including the momenta of the jet, the integrations over piz are fixed by residues approx- 
imately to Piz ~ E —pf±/ C^E). The contributions from the poles of the static potential 
are suppressed by a factor oc exp{—^\). In this model we neglect the small coUisional 
energy loss associated with target: dE^i/dx ~ Aa'^T'^ log E / /j, ~ 0.3 GeV/fm [1, 2]. 

Elastic scattering modifies the jet distribution in eq. (19) in this eikonal limit as 

x{DnJ{p-Q).rip-Q')) 

||p<»'(p-e)/n(<*V,|g)(||r(,',) 

x{27,f5\Q^-Y,q,^), (21) 

which includes the expected jet acoplanarity [22, 23] resulting from random walk but 
no recoil energy loss. 

The major simplification in the above derivation arises from the assumption that 
1^ ~ ~ 1/-^^ <^ /i <C A so that the main source of non-diagonal dependence arises 
through the T{qi± — q/j_) factors. This is what allows us to recover the above classical 
elastic parton cascade form in terms of the effective differential cross sections 
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The elastic cross section is then aei = {CjiC'iii) / D A){A.T:a1/ ^'^) which for a glue jet 
propagating in a glue plasma with = 0.5 GeV is typically (Tgi ~ 2 mb. (See Ref. [24] 
for a new Monte-Carlo implementation of such an elastic parton cascade.) 

We will neglect the influence of elastic multiple scattering on the inclusive distribu- 
tion of jets because is broadly distributed and the elastic energy loss is relatively 
small. In this case, p^^\p — Q) ~ p^^\p) can be factored out of eq. (21) leading to 

pir^p) ^ P^'\p) i'^eimr : (23) 

which is simply proportional to the free space spectrum times a geometrical factor (the 
Glauber multiple collision probability). 



3.2 Gluon Radiation with Quantum Cascading 

The amplitude to emit a gluon together with the jet including final state interactions 
with Us scattering centers is given by 

M^y^^ ^Mn. = J d'P J(P)e^^-o J n (-2ipO)V-(g;)e-^>^'r(z) j 

X MnAk,p;qi,---,qns) s\p - k-p + ^^Qi) ■ (24) 

This leads to the conditional double inclusive jet and gluon probability distribution 

d^k d^p 



d'vi:i^p^z\k,p) 



(27r)32c; (27r)32po ' 



ris 1 

pt\Kp) = (IM"^^ ® •Mn.r) - (n 77)Tr (IM"^^ ®-MnJ% , (25) 

which involves averaging over all initial target colors and summing over all final colors 
and the locations Xi of the target centers. 

After integrating over the energy-conserving delta functions we arrive at the mo- 
mentum space distribution density 

ptiik^p) = / fl {2p'fvmv*m{c,{^/DA)^ 



X 



]Je-(M>^' } J{k + p-Y,q^r{k+p-Y,q;) 



i=l 



X 



Ti(MnXk.P\qu---.qn^MiXk.P\qU---.qL)) ■ (26) 



The amplitude Adns is a sum over all time-ordered graphs involving Us scatterings in 
which one gluon is radiated. 
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For a slowly varying source, J{p + k — Q) ^ J{p + k) can be factored out of the 
integral in eq. (26) leading to 



P^lihp) ^ \J{p + k)\'j\[ 



\ V(27r)3 (27r)^ 



Da 



-T{qi± - qL) 



Us 



X 



i=l 



X Tr(^Mns{k,p;qu---,qnjMi^{k,p;ql,---,q^j') . 



(27) 



The integration over the can again be performed by contour integration leading to a 
sum over many terms with longitudinal phases set by energy differences as in [8, 9]. In 
the eikonal limit, a simple shortcut to that calculation is possible when the scattering of 
potentials localized in the longitudinal direction at Zi can be approximated by impulse 
interactions at times, ti ^ Zi/c [9]. In that case, time ordered perturbation theory can 
be used to express eq. (27) as 



i{k,p;qi±,---,qn,±)] 



X 



X 



Tr 



I \m=0 1=1 



ik,p;ql^, - ■ ■ ,qn,±) 



.m'=0/'=l 



(28) 



where the sums go over the 2"""'"^ — 1 different time ordered diagrams that we list in the 
next sections. The reduced radiation amplitudes, M.ns,m,i , for a given Ug are labeled by 
two integers, < m < and < I. Here m specifies the time interval tm < t < tm+i 
in which the gluon is radiated and I specifies which of the 2""""* final state interaction 
patterns with centers m + 1, • • • , occurs. 



3.3 Reduced Radiative Amplitudes 

It is convenient to associate with each amplitude diagram, M.ns,m,ii a-n t^s dimensional 
binary array, a, that decodes the label / for a given m as follows: 

a = (ai = 0, • • • , = 0, a^+i, • • • , a„J^ ^ / = (^ a,2^)/2"^+^ , (29) 

where cXi = for \ < i < m while (jj = or 1 for m + 1 < i < Ug depending on the 
integer part of //2'~"^~^. This is the base 2 representation of I x 2"^+-*^ — On^ - ■ ■ oi so 
that < Z < 2"^~"^ — 1. For a fixed m, Z = corresponds then to no final state gluon 
rescattering, while / = 2"^^™ — 1 corresponds to the gluon rescattering with all ng — m 
remaining scattering centers after tm- 
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For example, Fig. 2 shows Al5,i,io as the time ordered diagram with Ug — 5 
scattering centers, where ghion was emitted from the jet between ti and t2 and / = 10 
corresponding to a = (0,0, 1,0, 1) encodes the fact that the gluon received impulses 
<1±3, Q±5 at times, and t^. 

The binary representation is useful as a mask to extract the correct kinematic 
variables automatically. The transverse momentum of the gluon between interactions 
i and i + 1 is simply 

Ki = k±- Yl ^jQj-L = -^i-l + ^iQi-L ' (30) 

j=i+l 

— * — * 

with Kn^, = k±, the final gluon transverse momentum. For the example in Fig. 2, 
Kq^ Ki = K2 = k±- q3± - q5± and K3 = = k± - q5±. 



M 



5,1,10 
to 



t, t u 



i V 



^qi,ai ^q2,a2 "q3,a3 



q4,a4 



P 

k,c 



Fig. 2. Example of the reduced M5^i^io radiative amplitude for = 5. 



The complete amplitude matrix can be constructed in a few steps. First, the time 
integration over the emission times between tm and tm+i and the Ug impulses gives rise 
to the factor 

where the energy shift at the emission point is denoted 

K^^(^^E?i^ (32) 

The apparent singularity at cUmi = = due to the energy denominator is cancelled 

for finite t^+i — tm < oo. The singular tns+i = oo phase is always suppressed by an 
adiabatic damping factor (not shown). The only singular amplitude is in fact, the one 
corresponding to the factorization limit in which final state interactions factor out. 
Including the eikonal gluon emission vertex 2igsen{Km)p^ ~ Qs^i. ■ Km, the special case 
for m = ^s, I — reduces to 

Mn.,n.,o = -2ig/-^-^e''-^^c{an^ • • • ai) = e^(*--*°)^Ato (an, ■■■a,) , 

Us 
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where A^o is given by eq. (10). All the other amplitudes remain bounded < 
C{tm+i — tm)- On the other hand, in the formal Iq — > — oo limit, all the 2"= — 1 
singularities at uJmi = become exposed. 

The overall eikonal phase, ^mii acquired by the gluon through final state interac- 
tions is 

where Q.„dj = d^rni/dtj. 

The amplitude for color exchange leaving the final gluon with color c = 1, ■ • • , is 
specified by a particular color matrix, Tmic- The scattering of the jet with center i brings 
in a color matrix Oj = Ta^iR) in the R representation. Gluons interact with center i, on 
the other hand with Ta^^A) in the adjoint representation. It is numerically convenient, 
however, to express the complete color structure involving gluon final state interactions 
in terms of commutators in the jet {R dim) representation. Label sequentially in 
increasing order all rig = Y^k^k color matrices, ak where CTfe = 1, by ci, c^g. Label 
all remaining Ug — m — Ug matrices with k > m and ak — by bi, ■■• , bn^-m-ug- Then 

Tmlc = (bn,-m-ng ' ' ' bij ■ ■ ■ [c, C^J , " " " , Ci (a^ • • • Oi) . (34) 

In the example of Fig. 2, Ti_io,c = a4a2[[c, 05], aajai- Eq. (34) defines a straightforward 
algorithm for numerical computation of color factors, 

Cmlm'l' = E ( "^rnlcTmn'c ) ' 

c,ai,--,an 

needed as weights of the different interference terms in eq. (28) . 

Finally, the reduced radiation amplitudes with final state interactions are given by 

Mn,,m,i = 2igs ^2 I 2. - e' 2- 1 e'^'^Tmic ■ (35) 
3.4 Interference Form Factors 

If the "width". R^ ^ A ^ l//i, of the target transverse to the jet direction is large 
compared to its thickness, L, (i.e., R± ^ L > X ^ V/^)) then the off diagonal 
dependence in \qi± — qi±\ ~ 1/-R± in eq. (28) can be ignored in all factors except 
T{qi± — q/j_) (see Sec. 4.1 for more detailed discussion on this). In this case, the 
interference terms arising from the distribution of scattering times, ti, involve the 
following form factors: 

Fmlm'l' = { e'^*-' j^giWl'^mi _ ^tmU^ml^j f^^-itm' +l'^m' I' _ e-^*-"^m'i'^ ^ . (36) 

For 771,771' < 7is, the form factors insure that Fmlm'l' / i^mi^m'i') have no singularities 
except for m, m' = n^, and control the magnitude of the destructive LPM interference 
effects. The diagonal form factors are given by 

Fmlml = 2 - 2 ( COs{{tm+l - tm)uJml) ) , (37) 
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except that FnsOnsO = 1 since in^+i = 00 + ie. The extreme hmits of quantum cascading 
correspond to (1) the incoherent Classical Parton Cascade limit and (2) the coherent 
Factorization limit. In the first case, all non-diagonal form factors vanish, while in the 
second case all but the m = Ug diagonal form factor vanish: 

Fmim'V — '2Sm,m'Si,i' Classical Parton Cascade , 

Fmim'i' = Sm,ns^m',ns Coherent Factorization . (38) 

Simple analytic expressions can be obtained for the form factors if the ensemble 
average over the interaction times is taken as in eq. (20). The form factors involve a 
sum of four terms of the form 



1 



nr 

7 = 1 ^ 



iXjiuj -\- ■ ■ ■ + uj„ 



(39) 



where Tj = tj — tj-i, ojj = ujj + - ■ ■ + ujn,. Here Xj = (tj) is the average distance between 
the j — 1 and j^^ scattering centers. This ensemble is characterized by the following 
distribution oi t — tk — to: 




At— A 



icuXj 



(40) 



/4(i-VA,); A (^-1)! ■ 

The geometry of the plasma is particularly simple if A^ = A, with 

( (tk -to))^kX, { (tj - to)(tk - to) ) = k(j + 1)A' , j>k. 

For more realistic ensembles corresponding to an expanding target with centers 
distributed according to p{z,t). The form factors must be evaluated numerically via 

/ ,.E;;.fe-..M \ = -L d,„e-<»^'") /'" *„^, dt, 

\ I J to •'to J to 

x{t, t') = Id i aeip{z = i, i ) , 
Jt' 

= e-x(oo.*o)^"(^'^°) . (41) 

Tig. 

Here x{t, t') is the average number of interactions a jet moving with the speed of light 
through the expanding target suffers between t' and t with aei given by eq. (22). Note 
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that is the probabihty that only out of a large possible N :$> 1 target interactions 
occur in the medium. In this more general case, the interactions are distributed as 

atk 

lit) = x{t,to)/x{oo,to) . (42) 
For a static uniform slab of thickness L with density Pq9{L — z) 

^ E;:i \j = III. dtn^e'^-^^"^ ■■■ dtie''"^' , (43) 

It is instructive to compare eq. (43) with eq. (39). In the slab geometry setting to — 

kL , , k(j + l) 

Thus, if we set X = L/ (ug + 1) 

{h)^kX, {t^tk)^k{j + l)X'y^^ j>k. 

For rig^ 1 the geometry of both ensembles is quite similar if set A appropriately. Near 
the factorization limit, where YljLi Xujj <C 1 we can expand 

^.E.=ife-*o)-. \^l + ^X^JUJ,-rX'J2Il Kj + ^)^j^k , (44) 

j=\ j=i k=l 

where the only difference between eq. (39) and eq. (43) is that r = 1 in the first case 
and r = (n + l)/(n + 2) in the second case. The deviations from factorization can 
therefore be reasonably well approximated by our simpler ensemble. On the other 
hand, for large Xujj, the slab form factors oscillate due to the sharp boundary while 
the form factors in the eq. (39) ensemble decrease monotonically. The deviations from 
the classical parton cascade limit is therefore more sensitive to the detailed form of 
the geometry of the medium. We will continue to use below eq. (39) to reduce the 
complexity of the calculation. 

Since the final pattern of jet quenching depends on the actual distributions and cor- 
relations among scattering centers in general, one should evaluate eq. (36) numerically 
for any actual application to nucleus- nucleus collisions [12]. 

3.5 Relative Conditional Probabilities 

It is tempting to identify as in eqs. (15,16) 

Pel KP+k) 
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as the invariant soft gluon number distribution associated with Ug interactions. How- 
ever, because the radiated gluon can rescatter in addition to the incident jet, the 
conditional double inclusive density, p^^ad ^ cannot be factored into a product of the 
elastic distribution and the gluon number distribution as in the case of QED. The 
problem is that the production and gluon rescattering amplitudes interfere, and p'J'^} 
includes the process where the gluon is radiated before the first collision and the Ug 
recoil target partons are produced by elastic scattering with the produced gluon rather 
than the incident jet. Therefore, we interpret i?^"'-* as a relative conditional probabil- 
ity density to find a gluon in the final state with momentum k given a final observed 
jet with momentum p and given that Ug target recoil partons are observed (with any 
momentum) . 

Let Pn^ be the a priori probability that there are Ug scattering centers or target 
recoil partons. The final jet quenched spectrum is then given by 

PJ{P) = E {p^'\p) + [ — P^'^iP + R^;^\k)] , (46) 

where p*^°^ (p) is the initial hard jet distribution in eq. (12), and where the (wavef unction) 
renormalization factor 

Zn.^l-r j—Rf^\k) (47) 

insures that the integrated total number of jets is conserved. It in fact corresponds 
to the first order in ag term in the expansion of the inverse Sudakov factor for the 
probability that no radiation accompanies the jet. In this paper we only compute, 
R^^^\k), which hmits the discussion to single gluon emission: R^^'\k) = i?("«'"s=i) {k). 

Unfortunately, even in the rig — Q case the single gluon approximation is not accu- 
rate because 

Zo = l + ^W(£;o///), (48) 

TT 

can be significantly bigger than one, indicating that multigluon emission from the hard 
vertex is necessary. Summing the double logarithms leads to a modified Bessel function 
multiphcity distributions as discussed in [19]. The DLA leads to [19], 

Z, = l/Io{r{E)), 

r{E) = ^/8C^\og(E/|^) , 

r L(r) r(E) , , 

In the general case of multiple gluon emission, as required in practical applications, 
we must generalize eqs. (46,47) as 



pAp) = E (P ) + E / n ^ P^'^ iP + E iki,---XM , (50) 
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where the (wavefunction) renormahzation factor becomes 

^n. = 1 + E / n ^ • • • , 4j . (51) 

ng=l 1=1 ' 

The implementation of this approach requires Monte-Carlo methods that will be dis- 
cussed in the subsequent paper [12]. 

Prom the above formulas, we can infer that the exclusive Ug gluon number distri- 
butions associated with Ug scattering is 

rfiV("-"«) 1 , , ^ 

^1 • • -^n. ' ^3-. = • • • , kn,) ■ (52) 

a-^ki ■ ■ ■ d'^Kng '^ris 

4 Case of One Scattering Center 

We consider in this section the simplest case of gluon radiation in hard processes 
followed by one rescattering. 

4.1 Reduced Amplitudes and Transverse Profile Dependence 

The three graphs that generalize the Gunion-Bertsch analysis [6] to this case are shown 
in Fig. 3. 




Fig. 3. Three contributions to the soft gluon radiation amplitude Mj® Mi in case of one 
scattering. 

Following BDMPS [9], the three radiative amphtudes are given by 

A^i,o,o = 2i^,^^(e^*4_e^*o||)^^c^ (53) 
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x(e**^^^ - e'*°^^) [c,ai] , (54) 
A4i,i,o = 2iy.^^(-e^*4)cai . (55) 

Here ti is the impulse time when qij_ transverse momentum is imparted to the jet 
relative to its axis p. The color charge flow is modifled by an additional matrix ai, that 
appears in different orderings relative to the radiation vertex c in these amplitudes. A 
new infrared singularity k± — qi± — emerges due to the dipole radiation associated 
with the exchanged gluon. However, the phase factors cancel at that point and the 
amplitude remains finite proportional to ti — to- 

To keep the somewhat cumbersome notation compact we introduce the following 
reduced amplitudes 

ff_k± ^ _ (fe-fe -4 Q'iJ± 

p ' '-^(nJ2-im) (I. _ n. - n n- ' 

f^l. y'^ Hn Hlm)± 

= H — Ci , ^ {ili2---im){h32—in) — ^{hi2—jm) ~ ^{hh—jn) ' 

which we call "hard", "cascade", "Gunion-Bertsch" , and "Gunion-Bertsch cascade". 
The time integrals involve energy differences 

uq = k\l2uj = l/tf , uj{nn...i^) ^(k-Qi^-Qi^ Qim)±/'^^ , 



that control the formation physics. A quantity with two indices is defined to be the 
difference of two quantities with a single index (note that an array in parentheses counts 
as a single index). Thus, for example tio =ti — to, ujiq = uji — cjq. 

To make contact with the Gunion-Bertsch bremsstrahlung radiation distribution [6] 
and facilitate the physics interpretation we rearrange the sum of the three amplitudes 
as follows: 

Ml = -2i5,e**«'^"el • {Haic + Eie^^^o"^" [c, ai] + C^e-'^M'^'-'^o)) ai] } . (56) 

Including an adiabatic damping factor e~^*% we see that for — > oo the pure hard 
radiation formula for the = case is recovered. 

lim Mi = ai Mo ■ (57) 

ti— »+oo 

On the other hand, in the to — > — oo limit. 

Ml ^ -2ig,ei_ ■ 5ie**i'^o [c, ai] = Mgb , (58) 
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corresponding to the isolated bremsstrahlung amplitude. This leads to the character- 
istic GB radiation spectrum 



dyd^k^ ^7r^klik-q,)l' ^^^^ 

that has to be averaged over the transfered momentum. Here only the adjoint 
Casimir Ca = Nc enters and the result is independent of Cr. Note that in this limit, 
the k± — qi± is exposed and in [6] is regulated with an appropriate hadronic form 
factor. 

For the theoretical application to a parton penetrating through infinite matter, 
this singularity is regulated by the in- medium self energy of the gluon, 11^^ (A;). The 
transverse plasmon modes satisfy uj"^ = k'^ + nr('^, k) in an infinite plasma medium, 
with UT{LOpi{k), k) pa ^ and /x^ = g'^T'^{2Nc + Nf)/Q is the Debye screening 
mass in thermal pQCD. For A;^ ^ T on the other hand, u!pi{k) ^ kz-\-{k\-\-ii^ /2)/{2kz). 
Thus the infrared singularities are regulated by a scale ~ /i^. Near the light cone, the 
one loop approximation for the self energy is not likely to remain accurate though. 
In fact, gauge invariance requires that the same ii^ screening scale regulate both the 
potential singularity and the above infrared divergences [8, 9]. 

— * 

In contrast to the above, ^ —oo limit, the k± — qi± singularity is also regulated 
by canceling phases for tio finite. In this case of interest here, it is the maximum of /x^ 
and uj/tiQ that regulates that singularity. For large uo ^ i^BH = most important 

for energy loss, the latter scale in fact sets the characteristic transverse momentum of 
the gluon in the k± > fj, domain. 

Summing over final polarizations and colors and averaging over initial colors, the 
Us — 1 gluon distribution from eqs. (28,53-55) is 

x/ttI E M^,rnAk,P;qi±)Mlrn'Ak,p;ql^)]) . (60) 

\ \l,m,l,m',l' Jit 

The main complexity in the above expression arises from the mixing of different qi 
and ql components. However, if the transverse profile T{x^) of the target is approxi- 
mately a constant over dimensions 3> l//x, then the difference \qi — qi\ ~ 1/-R can be 
neglected in some of the terms. The only sensitive dependence of that difference arises 
through the phase factors in eq. (53-55). A systematic expansion in that difference is 
possible changing variables to Q — {qi -\- q{)/2 and q — qi — q{ and expanding in q. To 
leading order in this expansion 

p'r^Ak,p) = ''9l\Jip + k)\^J^,^\Vm'J^,Tiq)x 
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{Tr [{Haic + BQe'''°^° [c, Oi]) • (Hcai + BQe-'''°^° [a^, c])) 
+Tr (|CQ|2e-^*io('^i-'^i) [c,ai] [ai,c]) 

(61) 

+Tr ((^oic + SQe^*io'^o [c, Oi]) • (7^6+^*1° ^'^i-'^") [oi, c]) 

+Tr ((^aic + 5Qe-**io'^o [c, oi]) • (7^6-^*1° ('^'1"'^°) [ai, c]) )^ . 

Note that coi ^ {k - Q - q/2f/2u;, uj[ ^ {k - Q + q/2f/2u;, and therefore, the 
extra phase factor multiplying the diagonal cascade contribution has a phase 

bi ■ q = tio{k - Q) ■ q/uj , 

where we drop OitiQii^ /u) terms. The integration over (Pq transforms T[q) therefore 
into transverse spatial profile T{hi) evaluated at the impact parameter 

— * — * — * 

hi = tw{k - Q)/u , 

at which a gluon produced at the origin with transverse momentum k — Q would pass 
by the scattering center at zi = cti. This tends to suppress the cascade contribution 
as intuitively expected for kinematic configurations where the gluon "misses" the tar- 
get. On the other hand, the diagonal hard and Gunion-Bertsch contributions remain 
unsuppressed since the rescattering of the jet occurs at zero impact parameter in the 
high energy limit. For a broad target, the transverse profile is approximately a con- 
stant T(6i) pa T(0) ~ I/ttR"^. However, in general the cascade term is modulated by a 
relative target transverse profile factor 

T{h) ^ T{h)/m . 

The relative conditional probability for gluon radiation in this Ug — 1 case for a fixed 
momentum transfer Q± — > therefore reduces to 

Rf\k,qiA_) = CR^{H^ + R{Bl + T{h)C^i)-R{H -BicosihoUJo)) 

-R T(6i/2) • Ci cos{twUJw) - 2Ci ■ Bi cositwUi)) } . (62) 

Note that the transverse profile factor modulates the interference term between C and 
H + B at the average impact parameter, bi/2 of the jet and gluon at the scattering 
center. Note that for tio — > there is a perfect cancellation of all but the hard radiation 
term since bi oc tio and T(6 = 0) = 1 . In fact, this factorization limit is recovered 
more generally when tio/uj 0. 

The first three terms corresponds to the diagonal terms which are independent of 
the time between the hard production point and the rescattering point. The next three 
terms are oscillating interference terms. The color factor R — Ca/ Cr is 1 for a gluon 
jet and 9/4 for a quark jet. 

It is interesting to note that if the momentum transfer qi in the scattering vanishes, 
then Bi = 0, C = H, and cui = luq and the above distribution reduces to the zero 
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scattering case in the previous section. This is a manifestation of gauge invariance 
requiring a finite qi to generate a color dipole that can modify the hard radiation 
pattern. 

It is clear that the classical incoherent parton cascade limit corresponds to assuming 
that all the oscillating interference term average to zero. In that case, the final gluon 
distribution is a sum of the initial hard radiation plus an isolated Gunion-Bertsch 

bremsstrahlung plus a rcscattcring "cascade" contribution reflecting the modification 
of the final transverse momentum spectrum due to rescattering of a gluon produced at 
to- That latter term broadens the transverse distribution of the emitted gluons. 



4.2 Relative Conditional Probability for = 1 

The magnitude of the interference effects in eq. (62) is controlled by dimensionless 
variables 



Note that this averaged interference cosine is small when > 1. This is the kinematic 
region where the relevant formation time, Tj = 1/cjj = 2uj/{k± — qi_L)'^, is smaller than 
the mean free path. For a fixed energy uj this restricts the angle of induced emission 
radiation to exceed 9min = ^/Vi'^^)- Below 9min complete destructive interference 
cancels all but the hard production amplitude. 

The time averaged relative conditional distribution for fixed can be written as 

R^^\k,qi^) = i?f (64) 

where Rf^ is given by eq. (15) and the modification factor due to = 1 rescattering 
is given by 

/-i^ ^ f K,"^ KKl K.Ak,\ ^ f K,^ kAk,\ 1 

:f^^^ = ^ + R{^ + 2^-2-^] -R\ — - ^ 



-R 



\kkij i+em-Hi? \ 



kki / 1 + ^^(k — ki)'^ \ kki kl J 1 + ^^ref 



2 (2 
K ^ ^2 J K in KKl 



^2+^^ |i+e/^2+^i+e«i 

/kAk\( 2kI {k-kif k' \] 

\kkj[i+eiii i+eia-ni? i+e~^v i ' 

where Ak = k — k± ■ qi^/ jj?' = {k^ + kI — 9\)/2. The tilde denotes ki = Ki + e^, that 
results when a plasmon mass shift, uo'^i = e^/x^, is introduced in all propagators. 

Eq. (65) shows that for finite ^, the k± = qi± singularity of free space perturba- 
tion theory eq. (59) is automatically regulated by the interference terms. For finite ^ 
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therefore we can safely take the e = hmit and write 



KKi 



1 1 + ^2^2 1 + ^^kI 



However, in the formal classical parton cascade limit with ^ — > oo (a; ^ i^bh) the 
apparent singularity is exposed unless e 7^ 0: 

Ki \Ki K 



[H^ + R[Bl + Cf)) /H\ (67) 



To compare the full quantum parton cascade modification factor eq. (66) to the in- 
coherent classical parton cascade result, we regulated the singular classical limit with 
e = 1 

Qi±) = ^S\k±, qi±, oo)l^^ . (68) 

In the opposite ^ — > limit, on the other hand, we see that the factorization 
theorem, — > 1, is recovered. 

Another interesting limit is k±^ ^ qi± ~ fi. In that kinematic range, the induced GB 
bremsstrahlung can be ignored, {Bf pa 0) but the cascade term survives {Cf — > H^). 
Therefore, in the fj, <^ k± < cu kinematic regime 

^«(^x,gl^) 1 + R. (69) 

This enhancement of the high k± tail arises because rescattering of the gluon is en- 
hanced by a color factor it! = Ca/Cr relative to the jet. As we shall see, this is a 
special case of the general asymptotic k\ ^ ns^J^ form in the case of collisions: 

^("^)(^^,5l^,.-.^„,^) '-'^'^^ {l + Rr. (70) 

Although it is possible to integrate over the momentum transfers analytically, the 
result even in the case of one scattering center is not instructive. Therefore, we average 
^^^•^ over qi±^ numerically for different values of ^ — ojbh/'^'- 

Note that we keep the large s/4fi'^ ^ 6EjetT/Afi^ upper bound in the integration finite 
because the kinematic boundaries include regimes where E"^ > > k^> s. 

In Fig. 4a the three components of the classical parton cascade envelop (^ = 00) for 
an initial quark jet are shown for the case // = 0.5 GeV as a function of the logarithmic 
"angle" log^Q(/t). The Hard factorization term corresponds to the unit curve. The GB 
induced radiation component dominates at small k and the cascade contribution is only 
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important at large k. The bump at ~ 5 is the remnant of the logarithmic singularity 
exposed in the ,^ = 00 limit, but regulated here with a plasmon mass ujpi = /i = 0.5 GeV. 

Fig. 4b shows the ensemble averaged quantum modification factor, F^^\ with 
IJ, = 0.5 GeV, A = 1 fm/c ( ljbh ~ 0.63 GeV) for seven different gluon energies 
uj. These correspond to ^ = 0.013, 0.031, 0.063, 0.13, 0.31, 0.63, and 1.3 that cover 
an energy range 0.5 GeV < a; < 50 GeV . In each case, the curves are cut off at the 
kinematic bound 

which corresponds to right angle radiation relative to the jet axis. Here (A/i.)^/4 1.5 . 
The curves illustrate well how the total destructive interference (LPM effect) for < 1 
reduces F^^^ to unity. In particular, the induced Gunion-Bertsch radiation component 
in Fig. 4a is almost completely suppressed due to interference with the hard production 
amplitude. Only the cascade term modifies the spectrum considerably at large k±. The 
onset of significant medium modification can be seen to follow qualitatively the point 
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Fig. 4. (a) The normalized time-independent envelopes and separate contributions to the 
conditional probability distribution of gluons associated with a single rescattering = 1 
of quark jet (Ejg^ = 50GeV) is shown as a function of the logarithmic "angle", log^g k^j jjp' 
with ^ = 0.5 GeV. These curves correspond to the incoherent parton cascade limit, (b) 
For finite tiQ = 1 fm/c case destructive interference Hmits the corrections to the hard 
self-quenched distribution to higher angles as shown for different gluon energies. 

It is convenient to define in general the reduced modification factor f^""'\K,^) 
through 

F(-)(«:,0^{l + ((l + it:)--l) ^f^-«:) , (72) 
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which isolates the kinematic boundary and the asymptotic color scaling factor. Note 
that the incoherent classical parton cascade envelope is not reached for any reasonable 
values of X/j, due the finite kinematic boundaries. 



5 Case of Two Scattering Centers 

The relative conditional probability R^^^ {k ) can be obtained through the averaging of 
the square of the matrix element AI2 on the scattering potentials in the qi and ^ 
space. 

In Appendix A we listed the matrix elements connected to the seven planar dia- 
grams, contributing to the soft gluon radiation. The sum of the matrix elements can 
be written as 

M2 = -2i^,e^*o'^° X el • {Ha2aic + ^ie^*i°'^»a2 [c, Oi] 

+ 4e^*2°'^° [c, 02] ai + 52(i2)e'(*2°'^°-*2i'^2) [[c, 02] , ai] 
+ (7ie-'*^°^i°a2 [c, ai] + (726-**^""^° Oi [c, 02] 

+ (7(i2)e^(*^°'^°-*2i'^^-*i°'^(i2)) [[c, 02] , oi] } . (73) 

In the limit where all ti — tj = 0, ti > tj, i,j = 1,2 one can easily see that all the 
terms but the "hard" scattering cancel 

M2 = -2i^,e**o'^» ex ■ #02010 . (74) 

The GB to —00 limit and the connection to the BDMPS [9] result we discuss in 
Appendix B. 

After calculating the appropriate color factors in eqs. (92)-(93) we obtain 

Tr {M2MI) = CIDr x{h^ + R [bI + Bl + + C^) 

+R'^ (-62(12) + ^(12)) } + Interference terms . (75) 

(See eq. (95) for the full result.) Wc can isolate the color prefactor together with the 
leading parabolic dependence ~ R^g^ in the form 

Rf\k, gix) = Rf\k)r^'\k^, q,^, glx, 6, 6) ■ (76) 

We note that when the interference terms drop at asymptotic /c5_ S> 2/^^ the enhance- 
ment factor scales as in eq. (70) 

^(^)(^x,gix,g2x) '^-^'^ {l + Rf. (77) 

The longitudinal ensemble averaging can be performed in cq. (75) over the times 
^10 = ^1 — ^0 5 hi — t2 — ti as indicated in eqs. (20,39). For the more general case 
Ai 7^ A2 we find 

( cos(a;„tio + -,^21) ) = (1 + + • (78) 
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The integration over the momentum transfer due to the rescatterings from the static 
centers is most readily performed in terms of the dimensionless variables 



« = 4' ^^ = % e.-^^^^l^, ^=1,2, (79) 



leading to an averaged quantum correction factor 

Jo (01 + 1)2 {02 + 1)2 Jo 27r 27r ' ^ ^ 

On Fig. 5a we plot the contributions to J-'^i - the incoherent parton cascade limit 
where all the exposed singularities in the propagators have been regulated by the Debye 
screening mass /i. At small transverse momenta of the radiated gluon the quantum 
modification is dominated by the "Gunion-Bertsch terms". In the limit of large k\/ ix^ 
the only terms that survive apart from the hard scattering are the "cascade terms". 
This represents the expected behavior of R^^\k) in the Bcthc-Heitler limit. It is 
important to notice that simple picture of the process is only possible in the incoherent 
case where one easily decouples production from rescattering and the separate terms 
acquire physical interpretation. 

Fig. 5b shows JF^^) plotted for the same set of gluon energies cu as in Fig. 4b (and the 
same values of i = 1,2). The picture of classical parton cascading is modified by 
the presence of interference terms that control behavior of the probability distribution 
in the intermediate LPM region. Due to the destructive interference one observes a 
smooth transition from the "hard term" at small A;^ (i.e. the factorization limit) to 
the Bethe-Heitler limit. 




Fig. 5. (a) The structure of the regulated {/j, = 0.5 GeV) contributions to the ^ ^ oo 
hmit (quark jet of Ejg|. = 50 GeV, Us = 2). (b) Characteristic behavior of F^^) for 
^ = 0.013, 0.031, 0.063, 0.13, 0.31, 0.63 and 1.3 . 
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It is remarkable that the detailed angular profile of the probability distribution is very 
similar (apart from the explicit rescaling factor) to the one in Fig. 4b. 

6 Power Law Scaling in and Estimated Energy 
Loss 

The amplitudes in the case of Us = 3 are tabulated in Appendix C. The main novelty 
at this order is the appearance of an interesting "color wheel" diagram derived in 
Appendix D. The amplitude, A^a is rearranged in eq. (115) in terms of the hard, GB, 
and cascade amphtudes. The diagonal (classical parton cascade) terms in with 
color factors (117) involve (1) the usual — l/kj_ with a relative color weight 1, 
(2) three simple GB induced radiation terms and three one scattering gluon cascade 
terms each with color weight of R = Ca/Cr, (4) three double gluon cascade terms and 
four GB-cascade terms with a weight i?^, and (5) one triple gluon cascade term with a 
weight E?. At small k^, only the three simple GB (5^ — > H"^) terms contribute since 
the finite cascade terms are negligible relative to the divergent hard and GB terms. In 
this small /c^ region therefore the term is simply amplified by a factor 1 + 3i? as 
illustrated in Fig. 6a. For large angles, on the other hand, all the GB and GB-cascade 
terms vanish while the cascade terms — > H^. Therefore, in that region the hard 
radiation is amplified by a factor {l-\-FCf. In between, near A;^ ~ /x all terms contribute 
and a bump that is sensitive to the screening scale Fig. 6a summarizes therefore the 
general scaling with Ug of the classical parton cascade angular distribution. The power 
law enhancement (1 + i?)"^ at large angles is the most important nonlinear feature 
since the bump near k ~ 1 is killed by destructive interference shown on Fig. 6b. 




Fig. 6. (a) The classical parton cascade limit 00) angular enhancement factor from 

eq. (72) for a quark jet (Ej^^ = 50 GcV, ^ = 0.5 GcV) vs \og^Q{k\/ ^'^). (b) The reduced 
modification factor including kinematic cut-offs for = 1,2 shows an approximate 
Us independence. 
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In the general case, the small k± enhancement, 1 + UgR, expected in the classical 
cascade limit also is cancelled by destructive interference as seen in Figs. 4 and 5. 
Here we note a remarkable simple scaling between the Ug = 1,2 cases. In Fig. 6b 
the reduced modification factors ^) extracted from F^"'^'>{k,,^) in eq. (72) are 

shown. The reduced factors are nearly identical in spite of the much more complex 
functional form of the R^g^ than i? W . 

The approximate scaling from Ug — 1 to Ug — 2 can be summarized as follows: 

R^J'^\uj,k^)^Rf\k^){l + {{l + Rr-l) f^'\K,0} ■ (81) 

The generic analytic structure of the reduced modification factor can be roughly ap- 
proximated by f^^^K,^) ~ + K^C^), which shows that the classical parton 
cascade result is approached only for k > 1/^, i.e., k± > yfu^jx For higher the 
validity of this scaling will be tested in [12]. 

Here we extrapolate eq. (81) to = 3 to make a rough estimate of the energy loss 
spectrum dl/du and the fractional energy loss AE/E to L ~ 3 fm. We emphasize 
that these estimates are included in this paper only to illustrate qualitative features. 
For quantitative estimates, multi-gluon showers that require a Monte Carlo analysis 
must be considered as deferred to ref. [12]. The normalization factor Zn^ » 1 for even 
modest ~ 3 in eq. (47) signals a high probability of multi-gluon emission. This is 
well known in the Ug = emphasized in eqs. (48,49). For example, Zq > 2 from 

eq. (48) for E > 10 GcV while DLA gives, Zq ^ 3. Eq. (49) then leads to (ug) ^ 1.2 
for this energy, implying that two gluon final states must be taken into account. 




Fig. 7. (a) The gluon radiation intensity dl/dw for a 5 GeV quark jet {ug = 0,1,2,3) 
based on (81) is illustrated, (b) The fractional energy loss is shown for the "self-quenching" 
and medium modified cases. 

With the above caveats. Fig. 7 illustrates the qualitative features of the angular 
integrated dl/du! and AE/E for a quark jet of energy 5 GcV penetrating a thin 
{ug = 1,2,3) plasma with = 0.5 GeV, A = 1 fm, ag = 0.3. As seen from eq. (15), 
the intensity distribution for the self-quenching (n^ = 0) case increases linearly as a 
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function of log(a;//x) and the fractional energy loss due to self-quenching is proportional 

to \og{E/^). 

With increasing n^, the approximate scaling in eq. (81) leads to: 

^^^^^^^^ 

Without the Z factor, the gluon distribution would scale approximately exponentially 
with plasma thickness at high k± > fj./y/^ due to the multiple collisions of the gluon in 
the target. However, the normalization factor compensates this increase to a very large 
extent as seen in Fig. 7. In fact, at this single gluon level of approximation, the energy 
loss is surprisingly small, ~ 300 MeV/fm, and only grows approximately linearly with 
increasing thickness L = n^A in spite of the large modification of the gluon angular 
distribution. 



7 Summary 

We investigated the angular distribution of radiated gluons and the energy loss of hard 
jet in thin plasmas, where the size of the plasma L is comparable to the the mean 
free path, A. The small number of rescattering terms made it possible to calculate 
analytically the full matrix element within the eikonal approximation. This paper 
established a systematic method for computation of all relevant partonic matrix ele- 
ments for higher Ug which are needed as input to a numerical Monte-Carlo study of the 
actual hadronic observables associated with jet quenching to be presented in a subse- 
quent paper [12]. We developed an efficient algorithm to automate the cumbersome 
color algebra and keep track of the many canceling phases involved in the non-abelian 
Landau-Pomeranchuk-Migdal effect. Detailed discussion for the Ug — 1,2,3 cases were 
presented to clarify the how the competing physical processes interfere in different kine- 
matic regions. An approximate power law scaling of the angular broadening emerged 
with the increasing number of the rescattering centers. However, the rapidly growing 
wavefunction renormalization was shown to limit the growth of final energy loss. For 
quantitative comparison with data, multigluon emission will be shown in GLV H[12] 
to dominate the jet quenching pattern at the observable hadronic level. 
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A Case of Two Scattering Centers 



In the approximation of a high energy incident parton and soft gluon bremsstrahlung 
the seven planar diagrams in Fig. 8 can be treated in the eikonal approximation spinless 
hmit. 



k,c 



^2,0,0 ^^^"^ 

^to ./"tl t2 



T 
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>.to t tl 
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to 



tl 

n — 



i t I 

I I 

^ (X) ^ (X) 

qi,ai^ q2,a2^ 



M 



2,0,2 
to t tl 
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Fig. 8. The contributions to Mj(g) M2 soft gluon radiation amplitude in case of two 
scatterings. 

We denote by the time at which the particle is produced in the medium and 
by tj, i = 1,2 the times of the successive rescatterings. In the approximation when 
A ^ /i"^ "backward" scattering can be neglected and the process is approximately time 
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ordered [9]. We denote by k the momentum of the radiated gluon and by g^, i — 1,2 

are the momentum transfers. In QCD with each vertex we associate a generator of 
the SU{3) group. We denote them by c for the radiated gluon and Oj, i = 1,2 for 
the momentum transfer vertices. The polarization of the radiated gluon is determined 
by €±. 

The corresponding amplitudes can be obtained with the procedure described in 
eqs. (30,33-35): 

M2,o,o = 2i^.^^r2^(e'*4-e'*4)a2aic , (83) 

At I 



2 



[k - qi)i 



el • (fc - gi)±^rt,^iz(^ 
x(e'*i^^ -e**°^^)a2[c,ai] , (84) 

^ ., (fc-g2)i .. (fc-g2)i , ^ ^ 

x(e'*i^^ - e'*°^^)ai [0,02] , (85) 
2,0,3 = 2z^s— ^ T2^e''2 20. e'*i 2- x 



(fc — g]^— ^2)^ S - - 

x(e'*i 2;:^ -e^*° ^::r^)[[c,a2],ai] , (86) 

A<2,i,o = 2z^,^^-^(e**^&-e^*i^)a2cai , (87) 



(fe-'?2)i , 



X (e^*^ - e'*' ) [c, 02] ai , 



-* 7 

At2,2,o = 2i^,^^-^(-e^*^^)ca2ai . (89) 

After regrouping the sum of the above amplitudes five distinct color factors remain 
in eq. (73). We label them as follows: 

Ci = a2aic , C2 = a2 [c, ai] , C3 = [c, 02] ai , 

C4 = oi [c, 02] , C5 = [[c, 02] , oi] ■ (90) 

In the case of QCD the "color" and "kinematic" parts of a diagram factor. Diagram- 
matic techniques have been developed to treat the color part alone [25] . 

We denote by Cr the Casimir of the representation of the incident parton. For 
SU {N) following the standard normalization for the generators we have 

-1 

Cf = Ca = N. (91) 
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The trace of the color factors can be taken separately in the square of the matrix 

element, yielding 



Diagonal pcirt: 



Off-diagonal part: 



Diagrams with 0: 



TrCiCl = C%Dr , 

Tr C2CI = Tr C3CI = Tr C^Cl = CaCIDr , (92) 
TrCsCi = C'ICrDr . 

TVCiCl = -^CaCI^Dr , 

TrC^Cl = TrCiC] = -i {Cr - ^Ca) CaCrDr , 

TrC2Cl = TrC2Cl = -\C%CRDR, (93) 
Tr C3CI = -TrCiCl = {Cr - ^Ca) CaCrDr , 
TtCaCt = \C\Cr . 

TrCiC] =TrC2C] =0. (94) 



With R = Ca/Cr, the relative conditional probability is given by 

Rf\k) = {h' + R (Bl + Bl + Cl + Cl) + R' (4^2) + ^^2)) 

R -> -> R ^ -> 

- R ( — ) Bl ■ B2 COs(t2l'^o) - R ( — ) Cl ■ B2 COs(t2l'^0 + hoLOi) 

+ 2R^ (7(12) • -62(12) cos(tioa;(i2)) -R{h ■ Bi cos(iioa;o) 

-\- H ■ Cl cos(tioc<;io) — 2Ci ■ Bi cos{tiQUJi) 



-R{1- —){H ■ B2 COs{t200JQ) +H ■C2 cos(i2o<^2o) 

- 2C2 ■ B2 cos(t20'^2)) - R (— ) (Ci ■ ^2(12) cos(tioa;i - ^21(^20) 

+ Cl ■ (7(12) cos(iioa;i(i2) - t2i<^2o) + Bi ■ ^2(12) cos(i2i<^2o) 

+ Bl ■ (7(12) cos(tiocc;(i2) + ^21^^20)) - R"^ {C2 ■ ^2(12) COs(tioCU2) 

+ C2 ■ (7(12) COs(iioa;2(12)) - B2 ■ i?2(12) cos(i2i<^2) 

- B2 ■ (7(12) cos(iioa;(i2) + t2i<^2)) } , (95) 
where we have used the concise notation introduced in Sec. 4. 

B The to ^ —00 Limit in the Case of = 2 

Going back to eq. (73) it is easy to check the to ~oo limit and make contact with 
BDMPS [9]. In the phases we will assume small adiabatic switch-off factors ±is when 
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taking the corresponding limits. 

The hmit of to ^ ""OO corresponds to removing the interference with the jet source. 
In this hmit only the two Gunion-Bertsch terms and the Gunion-Bertsch cascade term 
survive. 

M2 = -2igs e"l ■ {5ie**^"«a2 [c, ai] + 52e**^"° [c, 02] oi 

+ 4(i2)e^(*^"°-*^^"^)[[c,a2],ai]} . (96) 

This is the two-scattering generalization of the Gunion-Bertsch one scattering cen- 
ter formula [6] and corresponds exactly to eq. (2.22) of BDMPS [9]. However, when 
squared there are several extra terms compared to the one effecfAve term retained in 
the approximations employed there. In particular, from eq. (95) 



lim 00—^ = / d^kx ^ = ^Ca ( / d^k^ \bI + Bl 

to^-co duj J dyd^k± \J ^ 

+ RB^^-^2) - (-^1 • B2 cos{t2iujo) + Bl ■ ^2(12) cos(i2ic^2o) 

-2 52-52(l2)COs(t2ia;2))} ) . (97) 

— * 

After taking the kinematic limits to infinity and performing the substitution kj_ = 
kj_ — q2± in the Bi ■ B2 contribution one arrives at a result that in the large Nc limit, 
i.e. R = 2, the weighted energy distribution of the gluons has the form: 

lim u — — uj- 



to^-oo dw duj 
+ 



S^Ca (/ d'^k^ (B2 ■ 4(12) (C0s(i2ic^2) - 1)) 

l^j d^k±_ [Bi ■ 52(12) (cos(t2ia;2o) - 1 



^j^{2)fact. jjy{2)BDMPS 
= OU hCU ' ^ , 

dOi) duj 

- ^Ca d^k^ (Si • ^2(12) (cos(i2ia;2o) - 1))) , (98) 

where direct comparison to eq. (2.31) in Ref. [9] shows that 

dN^BDMPS _ p CAas / f ,21: ( jk- q2)± (A:-gi-g2)± 
du - \J'"' [{k-q2)l {k-qi-q2f_ 

(C0s(t2l'^2) - 1)^ 
d:'k {B2 ■ 4(12)(cOs(t2lC^2) - 1)) ■ (99) 

If we also require that t2 — ti ^ 00 then the Gunion-Bertsch cascade term is also 
vanishes and we reproduce two independent Gunion-Bertsch scatterings with gluon 
radiation. 
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C Case of Three Scattering Centers 



The interaction of the incident high energy parton with three scattering centers accom- 
panied by emission of one gluon is described (in the eikonal approach) by the following 
fifteen diagrams: 
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Fig. 9. The contributions to Mj(g) M3 soft gluon radiation amplitude in case of three 
scatterings. 

The corresponding matrix elements are easy to obtain through the procedure described 
in Sec. 3.3 

MsAO = 2i^,^ip^(e^*ife-e^*°fe)a3a2aic, (100) 
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(fc-gi)i (fc-9i)i 

x(e**i 2u, -e'*° 2- )a302 [c,ai] , 
x(e'*i 2- -e**° 2u. )a3ai[c,a2] , 



2ig. 



e±- {k — qi — 02) ± ■■^^ ''i-<''-'?2)i .^^ (fc~g2)i-(fc-';i-g2)i 



{k-qi- ?2)i 

(fc-<31-g2)'i ., (fc-'?l-'32)^ 

x(e**i 2^ e**o ^ )a, [[c, 02] , ai] , 



x(e'*i 2U. -e'*o 20. )a2ai[c,a3] , 



(A; - gi - qsJi 
x(e'*i 2^ e'*o 21, )a2 [[c, as] , ai] , 



{k-q2- qs)! 



e_L ■ (A; - 52 - g^) 1 ^■^^ fci-(^-g3)i ^.^^ (fc-g3)i-(fc-g2-g3)l 
x(e'*i 2^ e^*o ^ )ai [[c, a,] , 02] , 

2^^. "^ " ^ ~ ^ ~ ^)^ c^^3 '' X 
(/c - gi - 52 - qs)! 

. (k-q3)']_-(k-g2-q3)\^ (fc-g2 -93)1 -('=-91-92 -93)1 

Xe 2 2u. 21^ X 

('s-91-92-93)^ (fc-91-92-93)^ 

x(e^*i e^*o ) [[[c, aa] , 02] , a^] , 



2igs ^^J^^ {e'^^ 2i _ e^*i 2i jaaaacai , 



e-l-(/c-g^2)± ,,, ''i-(!;^2)l 

. ('=-92)1 V ('=-92)1 

x(e"2 2a. -e"i 2a. )a3[c,a2]ai, 



el • (A; - g3)± ^it2^iz(±i23>l 

., ('=-93)1 V (fe-'?3)l 

x(e**2 2a. -e'*i 2a. )a2[c,a3]ai, 



[k - q3)i 
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{k-q2- qs)! 

x(e**2 5= _e^*i ^ ) [[c, as] , 02] ai , (111) 

M3,2,o = 2iy,%^(e^*^fe-e^*^fe)a3ca2ai, (112) 



2 

^ 2 „ x2 



x{e'''^^ - e'''^^) [c, 03] 0201 , (113) 



k 



e • k 

M3,3fi = 2ig, ^ ^ {-e''^^)ca3a2ai . (114) 

Organizing the sum of the matrix elements in groups with common color and phase 
factors we obtain 

M3 = -2igse^'°'^° X el • {#0302010 + Bie^''°'^°a3a2 [c, oi] 
+ 326"'^°^° as [c, 02] ai + 536^*^°'^° [c, 03] 0201 
+ S2(i2)e^(*^'"^°-*-'^^)a3 [[c,a2],ai] 
+ 4(i3)e*^*^°"°-*^^"^)a2 [[c,a3],ai] 
+ 4(23)6^^*^°"°-*^^"^^ [[c,a3],a2] ai 
+ #(23)(i23)e-'(*^^"^°+*^^"(^^)°+*^°"°) [[[c, aa] , 02] , ai] 
+ (7ie-^*i°'"i°a3a2 [c, Oi] + (72e-^*2°'"^°a3ai [c, 02] 
^(5^g-ito0Q2^^ [c, 03] 

+ (7(i2)e-^(*i°'^(i^)°+*^i'^^°)a3 [[c, 02] , ai] 

+ qi3)e-*(*i»'^(i3)o+*3iu.3o)Q^ [[c, 03] , ai] 

+ (7(23)e-'(*^°"(^^)°+*^^"^°)ai [[c, 03] , 02] 

+ C(i23)e-*(*^^"^°+*^i"(2^)«-*^°"(i23)o) [[[c, as] , 02] , ai] } . (115) 
Eleven distinct color factors are present for Ug — 3. We label them as follows: 
Ci = 0302010 , C2 = 03^2 [c, oi] , C3 = 03 [c, 02] oi , 

C4 = [c, 03] 0201 , C5 = 0301 [c, 02] , Ce = 0201 [c, 03] , 

(116) 

C7 = 03 [[c, 02] , ai] , Cs = 02 [[c, 03] , ai] , Cg = [[c, 03] 02] Oi , 
Cio = ai [[c, as] , 02] , Cn = [[[c, 03] 02] , oi] . 
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Computing the quantum amplitude requires knowing the color traces. These traces 
can be organized into a symmetric matrix: 

Diagonal terms: 

Tr CrCl = C%Dr , 

TVC2CI = Tr C3CI - TtC^CI = TvC^Cl - Tr CgC^ = CaC^Dr , 

TtCtC^ = TrCsCl = TrCgC^ = TrCioCjo = CIC^Dr , ^^^^^ 

Diagrams with 0: 

TrCiCj = TrCiCl = TrCiCl = Tr CaCjt = Tr CaC^J = Tr CaC^ = , 

TrCaC^ = TrCgC^ = TrC^Cl^ = TvC^Cl = TrCjClo = TrCgCJo = . ^^^^^ 

Traces involving the "Wheel" diagram (see Appendix D): 

TrCiCio = TrCaCji = TrC^Cl = TrCrCl 
TrCiCji = TrCaCjo = Tr CgC^ = Tr CgC^ 

Remaining off-diagonal terms: 

TVCiCl = -^CaC'rDr , 

TrCiCl = TrCiCt = -i {Cr - ^Ca) CaCIDr , 
TrCiC] = TrCiC^ - -i {Cr - ^CaY CaCrDr , 
TrCaC] = TrCaCj = TrCsC] = TrCgC^ = -\C\CIDr 
T1C2CI = TVC2CI = TrCsCJo = TVC4CI = = TrCsCJo = -\{Cr- ^Ca) CICrDr , 

= {Cr - \Ca) CaCIDr , 
TrCsCi = T^CiCl = -TrCgCl = -Tr CgC^ = \CICIDr , 
TrCsCl = TrC4C| = TrCyC^ = TvCrC^ = -^C^CrDr , 
TrCiCl = {Cr - ^Ca^CaCrDr , 

Tr CiCl = TvCiCl = -Tr CqCI = -TrCeCl = H^fl " ^^a) C\CrDr , 
TrCiCl = -TrCsCl = -TiCsCl = ^C^CrDr , 
TrCgCJo = {Cr - ^Ca) C\CrDr , 
TrCgCli = -TrCioCji = ^C^CrDr . 

(120) 



= W- \C\CrDr , 

= -W+\C\CrDr. ("^^ 



37 



D The Color "Wheel" Diagrams 

In the case of three scattering centers some nontrivial color factors appear for the first 
time in eq. (119). They correspond diagrammatically to the "wheel" diagrams shown 
in Fig. 10. 




,t» (Si crs c>, 

I ^1 Wa 1 'I 
a 1 b 



Fig. 10. The "wheel" diagram in the fundamental and the adjoint representations. 
For the fundamental representation 

{Ta.Tfj} = ^5ab + dabcTc, (121) 

It is useful to introduce the complex tensor 

habc df^bc ~l" ^fabc i ^aftc ^bca ^ca6 ; ^aab — ■ (-^22) 

If the physical case of SU{3), we find 

= "^"^ {TaTbTcT(i) fkalflbmfmcnfndk = (^-:^5ab5cd+ -^h^ 

1 33 

X ^ {^5ab5cd + 5ac5bd + ^5ad5bc " Qdacrdrbd) = y (123) 

for the fundamental representation and 

W^j4 {fkalflbmfmcnfndk){fraufudtftcsfsbr^ 

= ^ {5SabScd + SacSbd + 56adSbc " Qdacrdrbd) X 

X| {56adScb + ^ac^ferf + 5>6abSdc " Qdacrdrdb) = 189 . (124) 

for the adjoint representation. 
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9 
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33 
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135 
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405 
2 




24 


189 


840 


2775 



Table 1. The A^^ dependence of the wheel diagram to be found numerically. 
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